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Abstract 

In present note the arguments in favour of high statisics measurements of the form- 
factors that describe 7r-meson decays into a lepton-antilepton pair plus a photon are 
given. It is shown that these formfactors may contain an important information on 
the dynamics of quark motion inside a hadron. 
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1 Introduction. 

The aim of the present note is to discuss the arguments in favour of future high statistic 
measurements of 7r-mesons decay formfactors that contain the information on quark struc- 
ture of mesons. The decays of light 7r-meson such as it — > 7 + e + e~ and ^ —>■ 7 + e ± z/ (in 
what follows they would be denoted as ir° — > 7 + ee and it — > 7 + eu processes) are of a big 
interest from the view point of studying the dynamics of quark motion inside a hadron. 
The pioneering works on kinematics of these processes and the relation between the vec- 
tor Fy and the axial-vector Fa formfactors as well as on the relation with the width of 
7T° — > 77 process and also on the account of Inner Bresstrahlung amplitude conrtibution 
one can find in [|TJ. 

There is a plenty of different models based on current algebra, vector dominance model, 
chiral field theories and etc. (see, for instance, the reviews 0) that were proposed to de- 
scribe the formfactors dependence on a square of the invariant mass of a final state lepton 
pair. This variable can also be interpreted as a square of a 4-momentum Q 2 transfered 
from the hadronic block of the corresponding Feynman diagram to a leptonic pair. Only 
one particular prediction on a shape of a Q 2 - dependence of the decay formfactor, namely, 
on the appearance of a dip in a region of small values of x = Q 2 /M 2 variable ||, would 
be discussed here. It should be noted that the high statistic experiment done at Saclay 6 
years later had presented the data [|J (for more details see the reviews in || and |J) that 
may be interpreted as an experimental confirmation of the prediction done in || . Never- 
theless, the systematic errors quoted in |3| are rather high, so more precise measurements, 
especially in a region of small values of x, would be obviously very important. In this note 



it will be shown also that the analogous dip-effect may reveal itself in the formfactors of 
the electroweak decay 7r =l= — > 7 + e ± u. So, it would be very interesting to perform with 
good systematic errors new high statistic measurements of the behavior of F^o_> ie + e - (Q 2 ) 
formfactor in a region of small values of Q 2 as well as to collect the analogous data on 
F n ±_>^ e ±„(Q 2 ) formfactor. 

The prediction done in was obtained in the framework of the relativistic constituent 
quark model (see references in 0- [JTTJ] ) ) which make use of a covariant equation for two- 
body wave function [12j that was derived on the basis of 3-dimensional quasipotential 
approach to two-particle relativistic equations in quantum field theory (QFT) . 



The important point of this model is that its mathematical apparatus incorporates, 
in difference with ordinary QFT amplitudes, the bound state wave functions of quarks. 
Really, pion is a bound state of light quark and antiquark tighten together by forces 
caused by gluon exchange. Therefore, an application of Feynman diagrams (originally 
proposed in QFT for calculation of scattering amplitudes of particles that are free in an 
initial state) for describing the processes, that include the bound states in initial state, 
may serve, definitely, only as some perturbative model approximation applied in a region 
where the nonperturbative effects play an essential role. To this reason it is natural 
to expect that the consistent amplitudes to be used to describe the decay processes of 
mesons have to include the wave functions that take into account the bound state nature 
of mesons. 

Finally the prediction of the appearance of a dip is a sequence of three main features 
of the considered model: 

1. relativistic motion of quarks bounded in spin 7r-meson state; 

2. Standard Model (perturbative QCD and SM Feynman diagram technique) form 
of a quark propagator that enter the quarks interaction amplitude describing the photon 
and lepton-antilepton (e + e~ or e ± v) pair production in a final state; 

3. large value of a binding energy in a pion considered as qq bound state. 



2 Main formulas. 

To explain the above statements the main points of the analysis performed in U for 
7t° — > 7 + ee case would be sketched below. The schematic view ( not a Feynman dia- 
gram!) of 7r-meson decay processes that illustrates the corresponding invariant amplitudes 
M n ^f(P\qi, 92) at quark level (k\ and k 2 are quark and antiquark moments, while pi and 
Pi are the moments of the lepton and the antilepton respectively) is shown in Fig. 1. 

In Fig.l / denotes a final state, i.e. / = 7 + 11 for 7T°-decay and / = 7 + ev for 
7r ± -decay. These amplitudes may be parameterized through the decay formfactors in the 
following way (g x = p + + p_ = pi + pj): 



M n o^ ie+e -(q 1 ,q 2 \P) 
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Figure 1: The schematic view of 7r-meson decay processes 
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Here e u (q 2 ) is a polarization 4-vector of a real photon with the 4-momentum q 2 , ( e) is 
an electric charge of a lepton and Gp is Fermi coupling constant of weak interaction and 
V u d (~ cos(6 c )) is the element of Cabibbo-Kobayashi-Maskawa matrix. Also the following 
notations are used for two ("orthogonal' 1 



to each other) structure tensors and A^: 
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A^{ qi \P) = g ta ,(P-qi) -P^qu 
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that define , respectively, the vector and axial formfactors. Both electromagnetic (EM) 
and electroweak (EW) final state currents can be defined by one and the same formula: 



with the factor (V — A) defined as 

' V = 1:A = 



V -A 



for e + e , 
V = 1 ; A = 7 5 ; for e ± v . 



(5) 



(6) 



To make clear, why the same arguments that were used in may be valid also for 
the Fy(q\) and F^qf) formfactors of 7r ± — > 7 + ev decay, let us start with the definition 
of the amplitude of the process. 

An amplitude M n ^f(qi,q 2 \P) of a 7r-meson (considered as a qq bound state) decay 
process may be presented in a framework of the relativistic quark model as a convolution 
(with the relativistic invariant differential volume element of the momentum space 4^rr) 

of a covariant 2-body bound state (B) wave function ^/ CT1(T2 5p(k 1 ) and the final state 
interaction amplitude T q g_,f, taken in a form of Feynman matrix element (see Fig.l). 

Let us consider a general case of ir — > ^*V* decay, where by V* the virtual photon 7* 
or virtual W* boson (see Fig.l) is denoted. The transition amplitude is defined, following 
the guide line of [|ll[] and 0], like 
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1 r d 3 ~ki 

M 7T o^ rV *(q u q 2 \P) = J -^-^ ai(72BP ('k 1 )T^Z 2 r v*( ( l^ ( l2;k 1 \P) ■ (7) 

In (7) the summation over quark polarizations a±, a 2 is supposed. Now one may define 
the formfactor of ir° — > 7*7* transition as follows 

M w o_ +rr (q 1 ,q 2 \P) = F 1T o__ rr (ql,ql)e l ie v 2 e pLl/pa q p 1 P C7 = 

= F n o^ rr (ql q 2 2 )ete v 2 VMP) , (8) 



where and e 2 are the the polarization 4-vectors of two virtual photons with the 4- 
moments (off mass shell) q\ and q 2 respectively. 

After these definitions one can determine the formfactor of ir° — > / ~fe + e~ decay as (see 
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Fn°^ye+e-(Ql) = -^0^7*7* (Qi > 0) • (9) 

The wave function in (7) is a solution of a covariant two-body equation [[J, 0, |T2| 
that has a three dimensional form due to the use of a covariant single-time method of 
describing of a relative motion of quarks in a system where 7r-meson has the 4-momentum 
P. It should be mentioned that the relativistic wave function is connected with the vertex 
function r(P|k 1; k 2 ) according to formula 

(we take the masses of quark and antiquark to be equal to m q ). The vertex function 
is used in diagram technique to include the interactions between particles, in our case 
quarks, which lines in Feynman diagrams do enter this vertex. Formula (10) allows to 
establish a more close analogy with the Feynman diagram approach and with a quark 
triangle diagram often used to describe the decay process. 

Let us note that in perturbative leading order the vertex function, or the triangle 
Feynman matrix element (corresponding to decay amplitude), is set to be T(P|ki, k2)=l. 
It should be mentioned also that the three-dimensional nature of the integration in (7) 
over the 3-vector kx moment components (as well as a three-dimensial form of the wave 
function in (8)) is caused by a passing to a single-time formalism ||13|| , [14fl and by the 



fact that in difference with the ordinary Feynman diagram technique, where the virtual 
particle moments are "off the mass shell", in 3-dimensional approach the momenta of 
particles are on a mass shell, i.e. p 2 = m 2 , but the equations for the bound state wave 
function are written "off the energy shell" like in "old-fashioned" perturbation theory. 
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The interaction amplitude * v *(qi, k\\P) has the standard QFT form: 

T^ rV ,{q 1 ,q 2] k l \P) = 



47r v /ggv4 2 (fc 2 )ei(£;i - + m q )e 2 {V - A)iff(fci) 
(ki - qi) 2 - m 2 



(<&<-<Z 3 ). (11) 



Here e 2 = 7^e 2j ,(q , 2) with 62^(52), being the polarization 4- vector of a photon (that 
would be finally treated as a real one) with the 4-momentum q 2 , while e\ = 7 Q eiQ,(gi) 
with eia(q'i) being a polarization 4-vector of a virtual boson (photon or W^) with the 
4-momentum gi. The value s q = ^frT c ■ J2 e l includes a number of colors n c and the 
summation is done over the squared charges of quarks appearing in a fermion loop of 
a diagram shown in Fig.l. Factor a is the electromagnetic coupling constant while g is 
equal to s/a in a case of n° — > 7 + e + e~ decay and g = (eV u d) /2y/2sin(6 w ) in a case of 
the process n — > 7 + ez/. 

Factor (V — A) is defined by (6) and it takes into account the structure of the vertex 
(see Fig.l) corresponding to an intermdediate V* (= 7*/W /± ) boson coupling to quarks. 

The spin structure of ir- meson wave function is taken according to || , |TT | as follows: 



MWp(ki) = ^(fe 2 )7 5 <^i) 2 ^ ; (12) 

where <j\ and a 2 are quark polarizations and 0sp(ki) is taken to be a scalar function be- 
cause in what follows we shall consider qq s-state (i.e. with zero orbital angular momentum 
/ = 0). It should be mentioned that if one shall put (12) (with setting ^p^jm = ^) 
(7) and then substitute the wave function by the vertex function according to formula 
(10) and take there r(P|k 1; k 2 )=l, then an exact expresion of QFT Feynamn matrix 
element would appear under the sign of the integral. The nature of a 3-dimensial form 
of the integration can be easily understood on the basis of widely used rather straitfor- 
ward way of passing to a 3-dimensial formalism. In this approach one starts with the 
expression of a decay amplitude taken as an integral convolution ( with the 4-dimensional 
integration volume element) of a two-time 4-dimensional Bethe-Salpeter wave function 
with a 4-dimensional Feynman amlitude. Then by performing the subsequent equating of 
fermion and antifermion individual times in Bethe-Salpeter wave function by intoduction 
of (^-function, having the difference of these 2 times as its argument (see []13"1|- [fig] ) one 
gets a 3-dimensional equation in the momentum space. 



3 Dip-effect in n — ► 7 + e + e decay formfactor. 

In a case of tt° — > 7 + ee only the V=l term in (11) is taken in the amplitude, so it has a 
pure quantum electrodynamical QED form. 

After substituting of such a final state interaction amplitude into (7) one has to perform 
the summation over spin polarizations, what leads to the appearence of the corresponding 
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trace of 7 matrices, including those that were summed up with two polarization 4- vectors 
e iA»(<Zi) an d eivi^i) in (H)- The Lorentz indexes of these 7 M and 7^ matrices, associated 
with two fermion-boson interaction vertices shown in Fig.l, define the Lorentz index 
structure of the expresion for the calculated trace which is equal to V^ v (qi\P) in (3). 

Two polarization 4-vectors of photons ei M (gi) and e 2j ,(g 2 ) are not included, according 
to the definition (8), into the expression for the decay formfactor. The last one, thus, 
is a Lorentz scaler and is defined only through the quark block of the diagram shown in 
Fig.l, that includes only two quark wave function and quark variables which enter the 
amplitude (11). 

After calculation of the trace and the separation of its convolution with photon polar- 
ization vectors from the ampitude one may pass to performing the integration over the 
angular variables. Those are left, in a case of s-state, only in the denominator of quark 
propagator in (11) 



1 = _J L_ (13] 

{h-qi) 2 -m 2 q 2q-k 1 A + z' 1 ' 

where 

A-*-"**. (14) 

2q ■ fci V 1 

— * —* 

Here z = cos{9) = (q\ ■ k±)/qki and the notations k\ = \k\\ ; q — \q\\ are used. The 
following relations are valid also: for a real photon momentum we have q\ = q% — {(h) 2 = 0, 
i.e. q® = 1 <?2 1; for quark 4- momentum k\ = m 2 , because, as it was mentioned before, within 
the approach used in M the 4-moments of particles are "on the mass shell" but out of 



the covariantly defined "energy shell" [p~2 



We denote the square of the 4-momentum qi—(qi, q\) of a virtual vector boson V that 
produce a final state lepton-antilepton pair as q\ = Q 2 , keeping for the modulus \q\\ a 
notation q = \q\\. Thus, A = (Q 2 — 2qlk1)/(2q ■ k\) and we get an expression 



m q s q a I f dk\ ■ k\ (j)Bp(ki) f +1 dz 

^24. 
2q~~k x 



F^ ri ( qi ) = -=— 2tt / -— — — - ■ / } . (15) 



2vrM^ J 2kl 2q ■ k x J-i q 2 - 2qX , 

— cTZ — E r Z 



The appearance of the factor q = \q\\ in the denominator of (13) and finally in (15) 
has an important sequence that, possibly, may be experimentally observed. Really, due to 
the relation q\ = (P — qi) 2 , following from the 4-momentum conservation law P = qi + g 2 
we get a relation M 2 — 2Pq\ + q\ = 0. This invariant formula can be rewritten as 

2Pq x = M 2 (l +x), x = Q 2 /(M n ) 2 , (16) 

where from one can get in the pion rest frame (P = 0) the relations for the components 
of the 4- vector qi = (q®, q\) (keeping in mind our notation q = \q\\ and the definition 
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<f = \qi\ 2 = qf - (oi) 2 )- 




Ml + Q 2 
2M n 



M n (l + x)/2, 



(17) 



q = ml = 



Ml - Q 2 
2M„ 



M n {l - x)/2. 



Thus we see that after integration of the propagator in (15) over z = cosO one gets 



So, a possible growth of F(x) in the region of x ~ 1, if it can be observed in a data, 
may serve as a confirmation that the choice of the propagater in (11), as well as of the 
quark amplitude as a whole in a form of (11), i.e. in a standard for perturbative QFT 
form, may be quite a reasonable one and a consistent with data. 

A specific theoretical prediction of || , that is also connected with the the form of the 
propagator (13) from (11)( but is not comletely defined by it only) is about the formfactor 
behavior at small values of x . 

Firstly let us mention that the investigation performed in [jT7|j , where the so called 
"static" approximation for the wave function was used ( what is in fact equivalent to 
ignoring of the effects of quarks motion inside the meson) have shown that in this highly 
nonrelativistic approximation the slope of the formfactor has a positive value. 

In H another limiting case of ultrarelativistic quark motion inside pion (i.e. when to 
the integral (15) large values of k\ contribute mainly) was considerd. It was shown that 
in this limit the derivative of the formfactor has a negative value at sufficiently small x. 

Combining this observation, based on analytical calculations only, with the discussed 
above ~ (1 — x) _1 behaviour of the formfactor at x ~ 1, i.e. where it has a positive 
slope, one can suppose that for the relativistic bound state systems (and a light 7r-meson 
is a good candidate in such a case) described with the relativistic wave functions, the 
formfactor F(x) may have a minimum and, therefor, a changing sign of its slope. (It is 
worth mentioning that the present data for this slope, see, for instance, [[|, iclude the 
positive as well as negative values.) 

The origin of such possible prediction can be be understood from the analysis of the 
structure of formulae (15) without applying to any concrete form of the wave function. 

Really, the integration over ^-variable in (15) leads to an appearance of a logarithmic 
function |17|, (which may have different sign in different regions of its argument) under 
the integral sign 



the prediction |T7|,0 
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According to |§ the expression for the decay formfactor F n o^ +ee (x) being normalized 
to the constant of n° 27 decay 



F 7r o^ 7 * 7 (g 1 ) = f n o^ 2 -yF n o^ r * 1 (q 1 ) . (20) 



may be transformed to a form: 



1 Jo°° dXk4>(Xk) In \X(x, Xk)\ 
l-i T ' 4 So° d Xk<P(Xk)Xk 



= l-se-**(M p /m,-e-*») 

l_ xe x fc (M p /m g -ex fe ) ' ^ 

where 47i(j)(xk) = ^i0sp(ki) with ki = \ki\. The quark rapidity Xk — ln[{k® + k\)/m\ 
corresponds to the following parametrization of 4-vector components: 



k\ = m g chxk; h = |fci| = m q shx k - (23) 
The integral in the denominator of (21) defines the n° — > 77 decay constant |L7|, 



_ 32(2vr) 3/2 m (? as 
* 27 ~ M 2 



v 



L^2, = V tt, q q ■ / dx*0(Xk)x* • (24) 







It is clear from (22) that at x = one gets: -F(O) = 1 and thus 

F 7r o^ 7+ee (0) = 7^0^27- (25) 

Formula (21) presents the formfactor normalized to unity as a product of two factors. It 
may be also treated as a sum of two terms. One of them is (1 — a;)" 1 . It defines the 
monotonic growth of formfactor near x ~ 1 without a changing of the sign of a curve 
slope. 

The second one contains (besides the integral over the wave function multiplied by log 
of (22)) also the factor (1 — x) in the denominator and serves as a small correction to the 
main term (1 — x) -1 . In H it was shown that the sign of this correction depends on the 
sign of the logarithmic function ln\X\ under the integral sign and on what region of the 
integration of this ln\X\ over the rapidity dominats. The last circumstance depends on a 
shape of the wave function. 

The wave function <f>(xk) takes into account the bound state effect and, thus, has a 
nonperturbative nature. It serves in intgral as a weight factor for ln\X\ and defines what 
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region of quark rapidity may give the most contribution to the integral. Thus, the sign of 
the logarithmic function ln\X\ in this region would define the sign of a small additional 
integral term (appearing here as a correction to the leading term (1 — in (21)) and 
thus it defines finally the sign of the formfactor F(x) slope in a region of small values of 
x. 

In a case of M p /m q < 1 the numerator and the denominator in (22) are both positive, 
so the modulus sign in ln\X\ can be omitted. Then it is easy to check that for the values of 
quark rapidity Xk, satisfying the relation M n < 2m q ch(xk) , the numerator in (22) would 
be less than the denominator and thus the ln\X\ function would have a negative sign. So, 
we see that in a case when the binding energy of pion (as of the qq sytem )is negative ( 
i.e. it may be parametrized as follows: M n = 2m q cosj3) the condition < 2m q ch(xk) 
would be satisfied and thus the additional term to (1 — a;) -1 would be negative. It would 
lead, in principle, to a negative sign of the formfactor slope at small values of x. The 
value of a slope parameter depends on a shape of a particular wave function. 

From what was said above it is clear that the logarithmic function includes the infor- 
mation about perturbative amplitude (11). In this sense the log function in (21) fulfils a 
job of a perturbative probe by help of which one can get under the sign of the integral 
and test the shape of the qq bound state wave function in momentum space by means of 
variation the of external kinematic parameter x = Q 2 /M 2 value. 

Now after these general cosiderations it is a time to mention the results of where for 
numerical calculation of the formfactor behaviour two types of relativistic wave functions 
obtained in || as the solutions of 2-body relativistic three-dimensional equations with 
QCD inspired model potentials were used as well as one wave function that is an exact 
solution of relativistic oscillator model. All of these model wave fuctions have lead to 
a negative sign of the formfactor slope at small values of x (see below Fig. 2 where the 
result obtained within the relativistic oscillator model is presented). The results of other 
two QCD models give curves with the position of their minumums being from 20 to 40 
percents higher than that one of a curve shown in Fig. 2. 

In this connection it is worth mentioning the result of paper |18| where, for a sake of 
testing the method used in paper ||, a decay of purely QED system of muon-antimuon 
bound state into a real photon and e + e~ pair was considered. The invariant mass of two 
bound muons is much more higher than the value of two electron masses. So there is a 
more wide interval of Q 2 values may be attainable as comparing with a case of 7r-meson 
decay. 

At the same time the well studied apparatus of QED is known to be experimentally 
checked quite well and the tools for calculation of relativistic wave functions in a case 
of QED interactions in purely electromagnetic systems are also existing. So a muon- 
antimuon bound state is a good place to test the idea of the approach discussed above. 

The output of this kind of work done in |18| was finding out that the same dip in the 
formfactor of (/i/i) — > 7+e + e~ decay exists but its depth is of about one order smaller than 
that one found in F 7r o^ 7e + e - formfactor. This difference in a value of dip was obtained 
with the same expression for the amplitude (7) ( with A = ), but with the Coulomb-like 
relativistic wave function. Thus, from here one may make a conclusion that the depth of 
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Figure 2: The title of figure 2. 



the formfactor dip is defined by a value of binding energy of a system, which was found 
to be higher in a case of ir -meson within the theoretical models considered in ||. 

So, from comparison of the results of these two works, one may conclude that the 
value of the depth in dip, being measured at experiment, may give the information about 
the value of the binding energy in 7r-meson as a qq bound state system. 



4 Dip-effect in n —> 7 + e v decay formfactor. 

The detailed discussion of structure of the used expressions, perfomed in a previouse 
Section, allows now to pass easily to a case of n — ► 7 + ev processes. 

First let us note that the factors l)^r; 2)e; 3) j v {p +) pJ) that follow each other in the 

expression (1) do represent, respectively: 
1. the virtual photon propogator, 

g^/ql (26) 

It contains the metric tensor g a ^ which convolutes the Lorentz index of the 7°- matrix 
(that was included into t\ = 7 Q eiQ,(gi) in amplitude (11)) with the Lorentz index ji of 
lepton current jy(p + ,p_). This 7 - matrix was splited from the polarization 4-vector of a 
photon ei Q (gi) while extraction (according to formula (8)) of the formfactor F 7r o^ 7e + e - (qf) 
from the amplitude (7): 
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2. current j v (p + ,p^), that describs the lepton-antilepton pair production and is char- 
acterized by 

3. external to current j v (p + ,p_) factor e, i.e. an electric charge of electron - the 
QED coupling constant, not included into the expression for the current (5) for a sake of 
convention and for keeping the universal structure of current definition. 

All these three factors do correspond to the lines on the diagram of Fig.l that are 
external to the formfactor, which one, as it was discussed previously, contains only qq 
wave function and quark components of amplitude (11). 

Let us rewrite the formula (2) in an analogous way. For this aim we shall write the 
factors that have to appear according to the Feynman rules of Standard Model if we shall 
take them for the diagram shown at Fig.l: 

1.1 the virtual W-boson propogator 

2.1 current jv-A(P±iP») defined by (5) and (6), 

3.1 external to current 3v-a(p±iPv) factor of SM coupling constants, defined for a 
case of W exchange and the production of a pair of electron and neutrino in a final state 
/ = 7 + ev for 7r ± -decay as 

e/(2v / 2sm(^)), (28) 
while for the vertex of diagram in Fig.l, where W couples to quarks as 

{eV ud )/{2V2sin(9 w )). (29) 

The combination of all of this factors allows to write down the expression (2) in the 
following form 



M 7T ±^ e ± u (q u q 2 \P) 



MqD ■ V a M\P) + F A (qi) ■ A au {q 2 \P) 



x e 



ad 



W(^)), 



— • Jv-A\P±,Pv) e (92, 



(qi 2 ~ M, 



(30) 



w> 



which has the sructure analogous to formula (1) for QED process of 7r° — > 7 + e + e~ decay. 

If we shall consider in the last formula the limit q\ <C and take into account the 
relation 

(e/(2v / 2szn(^)) 2 = M^(G F /v / 2) 



(31) 

then we come to formula (2) that parametrizes the amplitude of tt^ — > 7 + e ± u process 
through the formfactors Fy and Fa- 
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Thus it is shown that the formula (30) has the same QFT structure as the formula 
(1), discussed in the previouse Section. Now if we shall write down the analog of formula 
(8) for a case of tt^ 1 —* 7 + e ± i/ decay and present the amplitude (7) as a sum of two terms 
according to two V and A parts of (V-A) factor in the amplitude (11), then defining by 
help of polarization vectors and e v 2 of virtual photon and the vitual IU-boson for each 
part the corresponding formfactor, we shall see that both of them would be defined by 
one and the same expression (21) as they are given by one and the same quark bloc like 
that one considered in a previuose Section. 

From here it is clear that the behaviour of the formfactors in a case of ir ± — ^ -7 + e ± t / 
decay must have the same dip as shown in Fig.2 for tt° — > 7 + e + e— process. 



References 

[1] V.G. Vaks, B.L. Joffe, Nuovo.Cim. 10 (1958)342; 

V.F. Muller, Zeit.Phys. 172 (1963)229; 173 (1963)938; 
S.G. Brown, S.A. Bludman, Phys. Rev.B 136 (1964)1160; 
P.De Baenst, J. Pestieau, Nuovo Cim.A 53 (1968)137; 
D.Yu. Bardin, S.M.Bilen'kii, Sov. J. Nucl. Phys. 16 (1973)311; 

D.A. Bryman, P. Depommier and C. Leroy, Phys. Rep. 88 (1982)152. 
L.G. Landsberg, Phys.Rep. 128 (1985)301. 

G. A. Kozlov, S.P. Kuleshov, V.V. Sanadze, V.I. Savrin, N.B. Skachkov, Zeit. Phys. 
C.21 (1983)63. 

H. Fonvieille et al, Phys.Lett.B 233 (1989)65. 

H.Fonvieille et al, Annales de Phys. CoU.n°2, Suppl. au n°3, vol.15 (1990)215. 

H.Fonvieille, These, Univ. Blaise Pascal-Clermont II, N° D'Ordre: E 420, PCCF T 
8908, 1989. 

N.B. Skachkov, I.L. Solovtsov, "Meson Form Factors and Covariant Three- 
Dimensional Formulation of Composite Model". Communication E2-11727, Dubna, 
JINR, 1978. 

N.B. Skachkov, V.I. Savrin, CERN-TH-2913, Jul 1980, Lett.Nuovo Cim. 29 
(1980)363. 

N.B. Skachkov, I.L. Solovtsov, Sov.Jour.Theor.Math.Phys. 54 (1983)116. 

V.N. Kapshai, A.D. Linkevich, N.B. Skachkov, V.I. Savrin, CERN-TH-3022, Jan. 
1981, Nuovo Cim.A. 66 (1981)45. 



12 



[11] V.I. Savrin, N.B. Skachkov, Sov. Jour.Theor.Mat.Phys.54 (1983)242. 

[12] R.N. Faustov, Ann. Phys.78 (1973)176. 

[13] E.E. Salpeter, Phys. Rev.87 (1952)328. 

[14] A.A. Logunov, A.N.Tavkhelidze, Nuovo Cim.29 (1963)380. 
V.G. Kadyshevsky, M.D. Mateev, Nuovo Cim.55 (1968)275. 

[15] R. Barbieri, R. Gatto and E. Remiddi, Phys. Lett.B 61 (1976)368. 

[16] L. Bergstrem, H. Snellman, G. Tengstrand, Phys. Lett.B 82 (1979)419. 

[17] L. Bergstrem, H. Snellman, Zeit. Phys.C 8 (1981)363. 

[18] G.A. Kozlov, S.P. Kuleshov, V.I. Savrin, V.V. Sanadze, N.B. Skachkov, 
Sov.Jour.Theor.Math.Phys. 60 (1985)647. 



13 



